Solar neutrino oscillations in the quasi-vacuum regime 
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Motivated by recent experimental data, we study solar neutrino oscillations in the range 
eVVMeV. In this range vacuum oscillations become increasingly affected by (solar and terrestrial) matter effects 
for increasing 5m?, smoothly reaching the MSW regime. A numerical study of matter effects in such "quasi- 
vacuum" regime is performed. The results are applied to the analysis of the recent solar neutrino phenomenology. 



1. Introduction 

Flavor oscillations either in vacuum Q or in 
matter Q provide a viable explanation to the 
solar neutrino problem The corresponding 
survival probability Py^^y^ is determined by the 
mass-mixing parameters bm? and uj, as well as 
by the electron density profile Ne{x) along the v 
path. 

For solar neutrino oscillations three character- 
istic lengths can be identified: the astronomical 
unit L — 1.496 x 10^ km, the oscillation wave- 
length in vacuum: 



in 



AnE „ 5m? 

= 2.48 X 10"^^^ km , 



Sm? 



E 



and the the refraction length in matter: 



in 



27r 



1.62 X 10^ 



km 



(1) 
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where Sw? is expressed in eV^, the neutrino 
energy E in MeV and the electron density 
iVe in mol/cm^. Usually, two asymptotic 
regimes can be identified: 1) the "Mykheyev- 
Smirnov-Wolfenstein" (MSW) oscillation regime 
H {5rr?/E > lO"'^ eVVMeV) characterized 
by L^sw _ 2.n,at < L, and the "vacuum" 
(VAC) oscillation regime g {Sm'^/E - 0(10"") 
eV^/MeV), characterized by L^^^ - L > L^at- 
In the MSW regime, the many oscillation cycles 
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in vacuum {L^f^ <C L) are responsible for com- 
plete decoherence of oscillations, and the survival 
probability depends on the detailed density pro- 
file in the Sun (and, during nighttime, in the 
Earth). Conversely, in the VAC regime the mat- 
ter effects suppress oscillations in the Sun and 
in the Earth, so that (coherent) flavor oscilla- 
tions take place only in vacuum, starting approx- 
imately from the Sun surface ||]. 

However, there is a transition regime [that 
might be called of "quasi-vacuum" (QV) oscil- 
lations], corresponding to the range 10~^° < 
5m? /E < 10"'' eV^/MeV, in which oscillations 
become increasingly affected by matter effects, 
and decreasingly decoherent, for increasing val- 
ues of Sm? /E. In this regime, characterized by 
Lmat < -^?sY — -^j none of the approximations 
used in the MSW or in the vacuum regimes can 
be apphed. In the past, QV oscillations have been 
considered only in a in few papers (see and ref- 
erence therein) since typical fits to solar v rates al- 
lowed only marginal solutions in the range where 
QV effects are relevant. However, more recent 
analyses appear to extend the former ranges of 
the VAC solutions upwards and of the MSW so- 
lutions downwards in Sm'^/E, making them even- 
tually merge in the QV range (see, e.g., (7|J|]). 

The QVO regime hase been recently studied 
in detail by Friedland in Q. Other numerical ||] 
and analytical ||l^ studies of QVO have also been 
performed. Here we discuss the results of |^ in 
the context of the most recent solar neutrino data. 
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2. Calculation of the survival probability 

From general quantum-mechanical arguments 
we can derive the probability that a v,, produced 
in the interior of the Sun is detected as a on 
the Earth 



P. 



0^© 



{1-Pq){\-P^ 



+ 2VPo(1 - Pq)P®{1 - Ps>) cos^ , (3) 



where Pq and are the transition probability 
Pv^^vi along the two partial paths inside the Sun 
(up to its surface) and inside the Earth (up to the 
detector), and the total phase ^ can be written as 



i = 



2E 



(4) 



where Rq is the solar radius and ^0 (^0) is the 
phase acquired in the Sun (Earth) matter. 

Equation ^ is a general formula that smoothly 
interpolates from the vacuum to MSW regime. In 
particular, in the vacuum limit Pq ~ cf, ~ Pq 
and ^0, ^0 — 0, we obtain the standard vac- 
uum formula. Conversely, in the MSW limit 
the oscillating term is averaged away and Pq ~ 
sin^ Lj^jPc -I- cos^wJ^(l — Pc), where uj^ is the 
mixing angle at the production point and Pc is 
the crossing probability P^m,^^^ (where is the 
heavier eigenstate in matter) , we obtain the stan- 
dard Parke's formula 

he evolution inside the Sun has been performed 
by direct numerical integration of the MSW equa- 
tion from solar the center to the solar surface.^ 
The electron density was taken from Fig- 
ure 2 in 1^ shows, in the mass-mixing plane, iso- 
lines of the difference c^—Pq (solid curves), which 
becomes zero in the vacuum oscillation limit of 
very small 5m? / E. The variable tan^ w has been 
chosen in order to chart the first two octants of 
the mixing angle range. In a semianalyt- 
ical approximation of Pq{'^ Pc) has been find 
by modifying the standard resonance prescription 
|14| with the Maximum Violation of Adiabaticity 
(MVA) prescription. This approximation gives 
Pq with a percent accuracy. 

^In all the QV range ~ tt/2 so that Pq ~ Pc. 
^In [pj has been explicitly verified that the results are in- 
dependent from the specific u production point, which can 
be effectively taken at the Sun center. 



The evolution inside the Earth is done by evolv- 
ing analytically the MSW equations at any given 
nadir angle ry, using the technique described in 
|l5|| , which is based on a five-step biquadratic ap- 
proximation of the density profile from the Pre- 
liminary Reference Earth Model (PREM) ^ 
and on a first-order perturbative expansion of the 
neutrino evolution operator. Figure 4 in |^] shows 
isolines of the quantity — Pq, which becomes 
zero in the vacuum oscillation limit of very small 
Svn? I E for two representative values of the nadir 
angle. 

Regarding the phases ^0 and ^0 in Eq. (^, it 
has been shown in Q that they are negligible for 
current fits. Only ^0 could have non-negligible 
effects in future high-statistics experiments. A 
semianalytic formula for the calculation of ^0 has 
been recently derived in | pO| . 

3. Decoherence of the oscillating term 

For increasing values of 5m? / E, the oscillating 
term in Eq. (H) become increasingly suppressed 
by energy smearing. The onset of the corre- 
sponding decoherence effect is best studied for 
almost monochromatic neutrinos, namely for the 
^Be and pep lines, having a width Ai? ~ O(keV) 
0. When Losc/L become greater of AE/{E) 
(where (E) is the average neutrino energy) the 
damping effect is important. This happen for 
dm^E > 10"* eV^/MeV (and before for con- 
tinuos spectra). 

For narrow line spectra, the energy-averaged 
oscillating factor can be written as (cos^) ~ 
D{5m^L/2{E)'^) cos^, where D is the modulus of 
the Fourier-transform of the spectrum |^]. Fig- 
ure 8 in ||l^] shows the function D as function of 
Sm^ /2{E) for the ''Be and pep lines. As antici- 
pated, damping factor is important for 5m^/E > 
2 X 10"* eVVMeV. For higher values of 5m'^/E, 
the oscillating term can be completely neglected. 

The survival probability P^^—^^^ has to be av- 
eraged in time, since it depends from i, both 
through Pq via 77 and through the phase ^ via 
the earth orbit eccentricity. While PQ{t) and 
cos^(i) can be be easily averaged separately 
|^,|l^, the calculation of the time average of the 
term Pq ( 1 — Pq ) cos ^ is in principle tedious. 
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Figure 1. 2i/ solutions to the solar neutrino prob- 
lem at 90, 95 and 99% C.L. of CL-f GA+K-KSK 
rates. 



Figure 2. solutions to the solar neutrino 

problem at 90, 95 and 99% C.L. of all data. 



Fortunately, the Earth effect is non negligible only 
when the oscillating term is doubly suppressed 
(Po ^ and £> - 0). Therefore, in Eq. (||) one 
can safely take ~ as far as the oscillatory 
term is non negligible. 

4. 21/ solutions to the solar neutrino prob- 
lem 

In this Section we apply the previous results to 
the analysis of the more recent solar neutrino phe- 
nomenology. Usually, the analysis is performed in 
the first octant in u) and splitting the 5m? range 
in the vacuum (few • 10^^^ < (5m^ < few • 10"^" 
eV^) and MSW {Sm^ > IQ-^ eV^) subranges. 

The lack of observation of a distorsion in the 
recoil electron spectrum in SuperKamiokande, to- 
gether with a rate observed in SuperKamiokande 
and Gallium experiments compatible with about 



1/2 suppression of the flux, strongly favor large 
mixing solutions. The LOW and VAC solutions 
tend to merge in the QV part of the plane, and ac- 
ceptable solutions appear also for w > 7r/4 
Consequently, the correct way to show the so- 
lutions is to use the entire mass-mixing plane 
without the artificial splitting in Sm? and for 
<Lu < 7r/2. 

Figure shows the combined analysis of the 
total rate information coming from the Su- 
perKamiokande (SK) Homestake |l9|, SAGE 
@], and GALLEX-GNO experiments. Fig- 
ure 1^ shows the combination between rate infor- 
mation and day and night spectral information 
from SK [Q. The technical details of the anal- 
ysis are explained in [2^ and references therein. 
The value of Xmin is 35.1 for the LMA, 37.8 for 
the LOW and 40.7 for the SMA solution, for 36 
degrees of freedom. 

The LMA and LOW solutions appear to be fa- 
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vored over the SMA solution in the global fit, al- 
though it is rather premature to think that the 
latter is ruled out. The vacuum solutions appear 
strongly disfavored by the nonobservation of SK 
spectral distorsions. From Fig. || we see that a 
solution to the solar neutrino problem emerges at 
99% C.L. in the QV range. The matter effects in 
QV solutions break the symmetry lo — + 7^/2 — u 
valid in the vacuum range (and that can be well 
observed in the vacuum solutions in Fig. |f|) . From 
Fig. ^ it appears that large mixing angle solutions 
are favored. For this reason, several experimental 
tests have been proposed to check them . 

5. Conclusions 

We have studied solar neutrino oscillations 
in the hybrid regime brn^ jE € [10^^°, 10"'''] 
eV^/MeV ( "quasi- vacuum" regime), in which 
the familiar approximations suitable in the vac- 
uum and MSW regime are not valid. We have 
discussed a general formula that interpolates 
smoothly between vacuum and MSW regimes. 
The effects of energy smearing on the oscillating 
term and the Earth matter effect have been taken 
in account. The results have been applied to the 
analysis of the most recent solar neutrino data, 
and suitable solutions to the solar v problem have 
been found in the QV range. 
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